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THE BOSE GAS: A SUBTLE MANY-BODY PROBLEM 

o: 

O ■ Elliott H. Lieb 

(N ; 

P-i- 

fl\ , Abstract. Now that the properties of the ground state of quantum- mechanical 

rf\ ' many-body systems (bosons) at low density, p, can be examined experimen- 

tally it is appropriate to revisit some of the formulas deduced by many authors 
4-5 decades ago. One of these is that the leading term in the energy/particle 
is 4irap where a is the scattering length of the 2-body potential. Owing to 
the delicate and peculiar nature of bosonic correlations (such as the strange 
A/ 7 / 5 law for charged bosons), four decades of research failed to establish this 
Oi , plausible formula rigorously. The only previous lower bound for the energy 

^J , was found by Dyson in 1957, but it was 14 times too small. The correct as- 

^_^ ■ ymptotic formula has recently been obtained jointly with J. Yngvason and 

Oi ' this work will be presented. The reason behind the mathematical difficulties 

(«_•? t will be emphasized. A different formula, postulated as late as 1971 by Schick, 

^_^ , holds in two-dimensions and this, too, will be shown to be correct. Another 

(__P ■ problem of great interest is the existence of Bose-Einstein condensation, and 

what little is known about this rigorously will also be discussed. With the aid 
of the methodology developed to prove the lower bound for the homogeneous 
gas, two other problems have been successfully addressed. One is the proof 
C^ ■ (with Yngvason and Seiringer) that the Gross-Pitaevskii equation correctly 

describes the ground state in the 'traps' actually used in the experiments. The 

S other is a very recent proof (with Solovej) that Foldy's 1961 theory of a high 

density gas of charged particles correctly describes its ground state energy. 



Or 



B 
> 

H " 1- INTRODUCTION 

c3 . 

Schrodinger's equation of 1926 denned a new mechanics whose Hamiltonian is 
based on classical mechanics, but whose consequences are sometimes non-intuitive 
from the classical point of view. One of the most extreme cases is the behavior of 
the ground (= lowest energy) state of a many-body system of particles. Since the 
ground state function ^(xi, ...,xn) is automatically symmetric in the coordinates 
Xj G B? of the N particles, we are dealing necessarily with 'bosons '. If we imposed 
the Pauli exclusion principle (antisymmetry) instead, apppropriate for electrons, 
the outcome would look much more natural and, oddly, more classical. Indeed, the 
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Pauli principle is essential for understanding the stability of the ordinary matter 
that surrounds us. 

Recent experiments have confirmed some of the bizarre properties of bosons 
close to their ground state, but the theoretical ideas go back to the 1940's - 1960's. 
The first sophisticated analysis of a gas or liquid of interacting bosons is due to 
Bogolubov in 1947. His approximate theory as amplified by others, is supposed to 
be exact in certain limiting cases, and some of those cases have now been verified 
rigorously (for the ground state energy) - 3 or 4 decades after they were proposed. 

The discussion will center around three main topics. 

1. The dilute, repulsive Bose gas at low density (2D and 3D) (with Jakob Yngvason). 

2. Repulsive bosons in a trap (as used in recent experiments) and the 'Gross-Pitaevskii 
equation (with Robert Seiringer and Jakob Yngvason). 

3. Foldy's jellium' model of charged particles in a neutralizing background (with Jan 
Philip Solovej) 

The discussion below of topic 1 is taken from Jit] (in 3-dimensions) and pj (in 
2-dimensions). That for topic 2 is taken from [§] (in 3-dimensions) and Q (in 
2-dimensions) . Topic 3 is from |5|. See also ||, 0]. 

Topic 1 (3-dimensions) was the starting point and contains essential ideas. 
It is explained here in some detail and is taken, with minor modifications (and 
corrections), from |6|. In terms of technical complexity, however, the third topic is 
the most involved and can be treated here only very briefly. 



2. THE DILUTE BOSE GAS IN 3D 

We consider the Hamiltonian for N bosons of mass m enclosed in a cubic 
box A of side length L and interacting by a spherically symmetric pair potential 

v{\xi-Xj\): 

N 

(2.1) H N = -fj,J2 A i+ J2 v (Wi-%\)- 

i=\ l<i<j<N 

Here x*i £ K 3 , i = 1, . . . , N are the positions of the particles, Aj the Laplacian with 
respect to Xi, and we have denoted ti 2 /2m by p, for short. (By choosing suitable 
units p could, of course, be eliminated, but we want to keep track of the dependence 



of the energy on Planck's constant and the mass.) The Hamiltonian (2.1) operates 
on symmetric wave functions in L 2 {K N , dx\ ■ ■ ■ dxyi) as is appropriate for bosons. 
The interaction potential will be assumed to be nonnegative and to decrease faster 
than 1/r 3 at infinity. 



We are interested in the ground state energy Eq(N,L) of (2.1) in the thermo- 
dynamic limit when N and L tend to infinity with the density p = N / L 3 fixed. The 
energy per particle in this limit 

(2.2) eo(p)= lim E (pL 3 ,L)/(pL 3 ). 

L — >oo 

Our results about eo(p) are based on estimates on Eq(N,L) for finite N and L, 
which are important, e.g., for the considerations of inhomogeneous systems in pi. 
To define Eq(N,L) precisely one must specify the boundary conditions. These 
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should not matter for the thermodynamic limit. To be on the safe side we use Neu- 
mann boundary conditions for the lower bound, and Dirichlet boundary conditions 
for the upper bound since these lead, respectively, to the lowest and the highest 
energies. 

For experiments with dilute gases the low density asymptotics of eo(p) is of 
importance. Low density means here that the mean interparticle distance, p -1 ' 3 
is much larger than the scattering length a of the potential, which is defined as 
follows. The zero energy scattering Schroedinger equation 

(2.3) -2^A-0 + w(r)-0 = 

has a solution of the form, asymptotically as r —> oo (or for all r > Rq if v(r) = 
for r > Rq), 

(2.4) i>(x) = 1 - a/\x\ 
This is the same as 

(2.5) a=limr-^4, 

where uq solves the zero energy scattering equation, 

(2.6) -2M«S(r) + v(r)u (r) = 



with Uo(0) = 0. (The factor 2 in ( |2.q ) comes from the reduced mass of the two 
particle problem.) 

An important special case is the hard core potential v(r) = oo if r < a and 
v(r) — otherwise. Then the scattering length a and the radius a are the same. 

Our main result is a rigorous proof of the formula 

(2.7) e (p) ~ infipa 
for pa 3 <C 1, more precisely of 

Theorem 2.1 (Low density limit of the ground state energy). 

(2.8) lim *^ = 1. 

pa 3 — >o Auppa 

This formula is independent of the boundary conditions used for the definition 
of e (p). 



The genesis of an understanding of eo (p) was the pioneeri ng w ork [ 11 ] of Bo- 
golu bov, and in the 50's and early 60's several derivations of ( |2.q ) were presented 
|l2| , [p.3| , even including higher order terms: 

(2-9) ^L = 1 + i^(pa 3 ) 1/2 + 8 (f - y/Z) (pa 3 ) log(pa 3 ) + 0(pa 3 ) 

These early developments are reviewed in [14]. They all rely on some special as- 
sumptions about the ground state that have never been proved, or on the selection 
of special terms from a perturbation series which likely diverges. The only rigor- 
ous estimates of this period were established by Dyson, who derived the following 
bounds in 1957 for a gas of hard spheres |l5| : 



(2.10) 1 <^L< l + 2Yl/3 



lOv^ ~ ^ppa ~ (1 - F x /3)2 
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with Y = Ait pa 3 /3. While the upper bound has the asymptotically correct form, 
the lower bound is off the mark by a factor of about 1/14. But for about 40 years 
this was the best lower bound available! 



Under the assumption that (2.8) is a correct asymptotic formula for the energy 
we see at once that understanding it physically, much less proving it, is not a simple 
matter. Initially, the problem prsents us with two lengths, a <C p -1 ' 3 at low density. 
However, ( |2.8| ) presents us with another length generated by the solution to the 
problem. This length is the de Broglie wavelength, or 'uncertainty principle' length 

4 = {pa)- 1 ' 2 . 

The reason for saying that i c is the de Broglie wavelength is that in the hard core 
case all the energy is kinetic (the hard core just imposes a ip = boundary condition 
whenever the sdistance between two particles is less than a). By the uncertainty 
principle, the kinetic energy is proportional to an inverse length squared, namely 
£ c . We then have the relation (since pa 3 is small) 

a < p~ 1/3 < t c 

which implies, physically, that it is impossible to localize the particles relative 
to each other (even though p is small) bosons in their ground state are therefore 
'smeared out' over distances large compared to the mean particle distance and their 
individuality is entirely lost. They cannot be localized with respect to each other 
without changing the kinetic energy enormously. 

Fermions, on the other hand, prefer to sit in 'private rooms', i.e., £ c is never 
bigger than p~ l > 3 by a fixed factor. In this respect the quantum nature of bosons 
is much m ore pronounced than for fermions. 

Since (|2.8| ) is a basic result about the Bose gas it is clearly important to derive 
it rigorously and in reasonable generality, in particular for more general cases than 
hard spheres. The question immediately arises for which interaction potentials one 
may expect it to be true. A notable fact is that it not true for all v with a > 0, 
since there are two body potentials with positive scattering length that allow many 
body bound states |l6|. Our proof, presented in the sequel, works for nonnegative 



v, but we conjecture that (2J3) holds if a > and v has no iV-body bound states 
for any N. The lower bound is, of course, the hardest part, but the upper bound 
is not altogether trivial either. 

Before we start with the estimates a simple computation and some heuristics 
may be helpful to make ( |2.8| ) plausible and motivate the formal proofs. 

With ua the scattering solution and fo(r) = uo(r)/r, partial integration gives 

|2 i „,|j- |2i i-: _ A _ I <r,.A„,' („\ („, („\ /_M2 i ../„M., I„\^2^ 



{2p\Vf \ 2 + v\f \ 2 }dx = 4tt / {2/*K(r) - (u (r)/r)} 2 + v(r)\u (r)] 2 }dr 

\x\<R JO 

(2.11) = 8Trpa\u (R)\ 2 /R 2 -* ^pa for R -* oo, 

if Mo is normalized so that fo(R) — ► 1 as R — > oo. Moreover, for positive interaction 
potentials the scattering solution minimizes the quadratic form in ( p.ll| ) for each R 
with Uo(0) = and uo(R) fixed as boundary conditions. Hence the energy Eq(2, L) 
of two particles in a large box, i.e., L 3> a, is approximately 8npa/L 3 . If the gas is 
sufficiently dilute it is not unreasonable to expect that the energy is essentially a 
sum of all such two particle contributions. Since there are N(N — l)/2 pairs, we are 



thus lead to Eo(N, L) rs 4irpaN(N—l)/L 3 , which gives (2.8) in the thermodynamic 
limit. 
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This simple heuristics is far from a rigorous proof, however, especially for the 
lower bound. In fact, it is rather remarkable that the same asymptotic formula holds 
both for 'soft' interaction potentials, where perturbation theory can be expected 
to be a good approximation, and potentials like hard spheres where this is not 
so. In the former case the ground state is approximately the constant function 
and the energy is mostly potential: According to perturbation theory Eq(N,L) sw 
N(N — 1)/(2L 3 ) J v(\x\)dx. In particular it is independent of p, i.e. of Planck's 
constant and mass. Since, however, J v(\x\)dx is the first Born approximation to 



87T/ua (note that a depends on /i!), this is not in conflict with (2.8). For 'hard' 
potentials on the other hand, the ground state is highly correlated, i.e., it is far 
from being a product of single particle states. The energy is here mostly kinetic, 
because the wave function is very small where the potential is large. These two 
quite different regimes, the potential energy dominated one and the kinetic energy 
dominated one, cannot be distinguished by the low density asymptotics of the 
energy. Whether they behave differently with respect to other phenomena, e.g., 
Bose-Einstein condensation, is not known at present. 



Bogolubov's analysis [11 presupposes the existence of Bose-Einstein conden- 



sation. Nevertheless, it is correct (for the energy) for the one-dimensional delta- 



function Bose gas |17[ , despite the fact that there is (presumably) no condensation 
in that case. It turns out that BE condensation is not really needed in order to 
understand the energy. As we shall see, 'global' condensation can be replaced by a 
'local' condensation on boxes whose size is independent of L. It is this crucial un- 
derstanding that enables us to prove Theorem 1.1 without having to decide about 
BE condensation. 

An important idea of Dyson was to transform the hard sphere potential into a 
soft potential at the cost of sacrificing the kinetic energy, i.e., effectively to move 
from one regime to the other. We shall make use of this idea in our proof of the 
lower bound below. But first we discuss the simpler upper bound, which relies on 



other ideas from Dyson's beautiful paper 15 



2.1. UPPER BOUND. The following generalization of Dyson's upper bound 
holds [§, jfij: 



Theorem 2.2 (Upper bound). Define p x = (N~l)/L 3 andb = (47rpi/3)" 1/3 . 
For nonnegative potentials v, and b > a the ground state energy of ( \2.j\ ) with peri- 
odic boundary conditions satisfies 

1_ « + (ii) 2 + I (2) 3 
(2.12) E {N, L)/N < Anrnpia b - u; 8 2 u; ■ 

I 1 — b ) 

For Dirichlet boundary conditions the estimate holds with (const. )/ L 2 added to the 
right side. Thus in the thermodynamic limit and for all boundary conditions 

(2rt\ & M l-YV*+Y*t*-\Y 

K ' ' Anppa ~ (1 - yV3)8 

provided Y = 4irpa 3 /3 < 1. 
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Remark. The bound (2.12) holds for potentials with infinite range, provided 
b > a. For potentials of finite range Rq it can be improved for b > Rq to 

1 _ (a) 2 + 1 (*) 3 
(2.14) E (N,L)/N<iTTfi Pl a- KbJ 2Kb ' 



;i-t) 4 



Proof. We first remark that the expectation value of fl2.l| ) with any trial wave 
function gives an upper bound to the bosonic ground state energy, even if the 
trial function is not symmetric under permutations of the variables. The reason 



is that an absolute ground state of the elliptic differential operator (2.1) (i.e., a 
ground state without symmetry requirement) is a nonnegative function which can 
be be symmetrized without changing the energy because (pj) is symmetric under 
permutations. In other words, the absolute ground state energy is the same as the 
bosonic ground state energy. 



Following [15 we choose a trial function of the following form 

(2.15) *(a;i,.. .,xn) = Ft(xi) ■ F 2 (xi,x 2 )- ■ -F N (xi,. .. ,x N ). 

More specifically, F\ = 1 and Fi depends only on the distance of xi to its nearest 
neighbor among the the points x\, . . . ,Xi-\ (taking the periodic boundary into 
account): 

(2.16) Fi(xi,...,Xi) = f(ti), U = min(|xi - Xj\, j = 1, . ..,«- 1) , 
with a function / satisfying 

(2.17) 0</<l, /'>0. 



The intuition behind the ansatz (2.15) is that the particles are inserted into the 
system one at the time, taking into account the particles previously inserted. While 
such a wave function cannot reproduce all correlations present in the true ground 
state, it turns out to capture the leading term in the energy for dilute gases. 
The form (p. 16) is computationally easier to handle than an ansatz of the type 



lli ■ f(\xi — Xj\), which might appear more natural in view of the heuristic re- 
marks at the end of the last subection. 
The function / is chosen to be 

(2.18) f(r) = [M/Mb) for < r < b, 

1 for r > b, 



with fo{r) — uo(r)/r. The estimate s ([2.12| ) and (2.14) are obtained by somewhat 



timate s ([ET 
as in |15|, 1 



lengthy computations similar as in |l5[, but making use of (2.11). For details we 



refer to W and 

A test wave function with Dirichlet boundary condition may be obtained by 



localizing the wave function (2.1 5|) on the length scale L. The energy cost per 



particle for this is (const. )/L ■ □ 

2.2. LOWER BOUND. In the begin ning it was explained why the lower 



bound for the bosonic ground state energy of (2.1) is not easy to obtain. The three 

plaj 

• The mean particle distance p~ x ' 3 . 



different length scales for bosons will play a role in the proof below 
• The scattering length a. 



• The 'uncertainty principle length' £ c , defined by fi£ c 2 = eo(p), i.e., £ c 



(pa) 



-1/2 
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Our lower bound for eo(p) is as follows. 

Theorem 2.3 (Lower bound in the thermodynamic limit). For a posi- 
tive potential v with finite range and Y small enough 

(2.19) ^>(1-CT 1 ' 17 ) 

with C a constant. If v does not have finite range, but decreases at least as fast as 
l/r 3+£ at infinity with some e > 0, then an analogous bound to (2.1f\) holds, but 



with C replaced by another constant and 1/1 7 by another exponent, both of which 
may depend on e. 



It should be noted right away that the error term —CY 1 / 17 in (2.19) is of 
no fundamental significance and is not believed to reflect the true state of affairs. 
Presumably, it does not even have the right sign. We mention in passing that C 



can be taken to be 8.9 | 10 | 



As mentioned in the Introduction a lower bound on Eq(N, L) for finite N and 
L is of importance for applications to inhomogeneous gases, and in fact we derive 



( 2. IS ) from such a bound. We state it in the following way: 

Theorem 2.4 (Lower bound in a finite box). For a positive potential v with 



finite range there is a S > such that the the ground state energy of (2.1) with Neu- 
mann conditions satisfies 

(2.20) E (N, L)/N > Airppa (l-C Y 1/17 ) 

for all N and L with Y < S and L/a > C'Y~ 6 < 17 . Here C and C are constants, 
independent of N and L. (Note that the condition on L/a requires in particular 



that N must be large enough, N > (const. )Y~ l < .) As in Theorem 2.5 such a 
bound, but possibly with other constants and another exponent for Y , holds also for 
potentials v of infinite range decreasing faster than 1/r 3 at infinity. 



The first step in the proof of Theorem 2.4 is a generalization of a lemma of 
Dyson, which allows us to replace v by a 'soft' potential, at the cost of sacrificing 
kinetic energy and increasing the effective range. 

Lemma 2.5. Let v(r) > with finite range Rq. Let U(r) > be any function 
satisfying J U(r)r 2 dr < 1 and U(r) = for r < Rq. Let 6cR be star shaped 
with respect to (e.g. convex with € B). Then for all differ entiable functions if) 

(2.21) I [p\V^\ 2 + iz#| 2 ] > pa f t/|Vf. 

Jb Jb 



Proof. Actually, ( 2.21 ) holds with p\W(f>(x)\ 2 replaced by the (s malle r) radial 
kinetic energy, p\d(j>(x) / dr\ 2 , and it suffices to prove the analog of (2.21) for the 



integral along each radial line with fixed angular variables. Along such a line we 
write (f>(x) — u(r)/r with u(0) = 0. We consider first the special case when when 
U is a delta-function at some radius R> Ro, i.e., 

(2.22) u(r) = J-. S (r-R). 
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For such U the analog of (2.21) along the radial line is 

(2.23) 



{H[u'{r) - {u{r)/r)f + ±v{r)\u(r)] 2 }dr > 







if i?i < R 



lia\u{R)\ 2 /R 2 ifi?<i?i 



where R\ is the length of the radial line segment in B. The case R\ < R is trivial, 
because /j,\dtp / dr\ 2 + ^v\ip\ 2 > 0. (Note that positivity of v is used here.) If R < R\ 
we consider the integral on the the left side of ( 2.23 ) from to R instead of R\ 
and minimize it under the boundary condition that w(0) = and u(R) is a fixed 
constant. Since everything is homogeneous in u we may normalize this value to 
u(R) = R — a. This minimization problem leads to the zero energy scattering 
equation ( |2.6| ). Since v is positive, the solution is a true minimum and not just a 
stationary point. 

Because v(r) = for r > Rq the solution, uq, satisfies uo(r) = r — a for r > Rq- 
By partial integration, 

(2.24) 
?. 
{n[u' (r) - (u (r)/r)} 2 + \v{r)\u {r)] 2 }dr = na\R - a\/R > fia\R - a\ 2 /R 2 . 



But \R — a\ 2 jR 2 is precisely the right side of (2.23) if u satisfies the normalization 
condition. 

This derivation of ( |2.21 ) for the special case (2.22) implies the general case, 
because every U can be written as a superposition of ^-functions, U(r) = J R~ 2 S(r— 
R) U{R)R 2 dR, and / U{R)R 2 dR < 1 by assumption. □ 

By dividing A for given points x± , . . . , xn into Voronoi cells Bi that contain 
all points closer to xi than to Xj with j =/= i (these cells are star shaped w.r.t. Xi, 



indeed convex), the following corollary of Lemma 2.5 can be derived in the same 
way as the corresponding Eq. (28) in [15]. 



COROLLARY 2.6. For any U as in Lemma [2J3 
(2.25) H N > fiaW 

with 



iV 



(2.26) 



W(x u 



,x N 



)=$>(** 



where tj is the distance of Si to its nearest neighbor among the other points Xj, 
j = 1, ...,N, i.e., 



(2.27) 



U(xi, . . . ,xjv) 



: mm \Xi 



(Note that ti has here a slightly different meaning than in ( 2.16 ), where it denoted 
the distance to the near est neighbor among the Xj with j < i — 1.) 

Dyson considers in |15| a one parameter family of U's that is essentially the 
same as the following choice, which is convenient for the present purpose: 



(2.28) 



U R {r) = 



3(i? 3 - Rl)- 1 for Rq < r < R 
otherwise. 
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We denote the corresponding interaction ( 2.26| ) by Wr. For the hard core gas one 
obtains 

(2.29) E{N,L)> sup inf /j,aW R (xi, . . . ,x N ) 

R (X1,..,I W ) 

where the infimum is over (xi, . . . , Xn) € A^ with \xt — Xj\ > Rq = a, because of 
the hard core. At fixed R simple geometry gives 

(2-30) irf W R ($ u ...,x* N )>(^--^) 

(xi,...,xjv) \R J PR J 

with certain constants A and B. An evaluation of these constants gives Dyson's 
bound 

(2.31) E(N,L)/N>^L=4*npa. 

The main reason this method does not give a better bound is that R must be 
chosen quite big, namely of the order of the mean particle distance p -1 ' 3 , in order 
to guarantee that the spheres of radius R around the N points overlap. Otherwise 
the infimum of Wr will be zero. But large R means that Wr is small. It should also 
be noted that this method does not work for potentials other than hard spheres: If 



\xi — Xj\ is allowed to be less than Rq, then the right side of ( 2.29 ) is zero because 
U(r) = for r < R . 

For these reasons we take another route. We still use Lemma 2.21 to g et int o 
the soft potential regime, but we do not sacrifice all the kinetic energy as in (2.25). 
Instead we write, for e > 



(2.32) H N = sH N + (1 - e)H N > eT N + (1 - e)H 



N 



with Tff = — J2i A» and use (2.25) only for the part (1 — e)H^. This gives 

(2.33) H N > eT N + (1 - e)^aW R . 

We consider the operator on the right side from the viewpoint of first order pertur- 
bation theory, with eTn as the unperturbed part, denoted Ho. 

The ground state of Hg in a box of side length L is ^(^l , . . . , xn) = L~ 3N / 2 
and we denote expectation values in this state by (-)o. A computation, cf. Eq. (21) 
in 0, gives 



4^(1-^) > (W R ) /N 



(2.34) > 4np (1 - i) (1 - ^) A (1 + 4tt P (1 - ±)(R* - i? 3 )/3)) . 

The rationale behind the various factors is as follows: (1 — -k) comes from the 
fact that the number of pairs is N(N - l)/2 and not N 2 /2, (1 - 2R/L) 3 takes 
into account the fact that the particles do not interact beyond the boundary of 
A, and the last factor measures the probability to find another particle within the 
interaction range of the potential Ur for a given particle. 

The first order result ( 2.34| ) looks at first sight quite promising, for if we let 



L — * oo, N — v oo with p — N/L 3 fixed, and subsequently take R — > oo, then 
( Wr) q/N converges to 47rp, which is just what is desired. But the first order result 
( [2.34| ) is not a rigorous bound on Eq(N, L), we need error estimates, and these will 
depend on e, R and L. 
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We now recall Temple 's inequality 18 for the expectations values of an operator 



H = Hq+V in the ground state (-)o of Hq. It is a simple consequence of the operator 
inequality 

(2.35) {H-E o )(H-E 1 )>0 
for the two lowest eigenvalues, Eq < E\, of H and reads 

(2.36) ^>(g)o- ( ^° 

provided E\ — (-ff)o > 0. Furthermore, if V > we may use E\ > E\ — second 

?(0) 



lowest eige nvalu e of H an d replace E\ in ( |2.36| ) by E\ 



From ( 2.34 ) and ( 2.36 ) we get the estimate 



(2.37) Eoi ^ L) >A^ap(l~£(p,L,R,e)) 

with 

l-£(p,L,R,e) = (l- £ )(l--i T )(l-M) 3 ( 1+ 4 fp(1 _^ )(jR 3_ ii ,3 ))) - 

(2.38) x (l- ^»- W ° 

V (Wr) (e[°> - im{Wr)o 



To evaluate this further one may use the estimates (2.34) and the bound 




(2-39) (W^o<3 R3 N R3 (W R )o 

which follows from U R = 3(i? 3 — R 3 ) )~ 1 Ur together with the Cauchy-Schwarz in- 
equality. A glance at the form of the error term reveals, however, that it is not 
possible here to take the thermodynamic limit L — > oo with p hxed: We have 
E\ = ett/j,/ L 2 (this is the kinetic energy of a single particle in the first excited 
state in the box), and the factor E{ — pa(Wfi)o in the denominator in ( J2.38J ) is, up 
to unimportant constants and lower order terms, ~ (eL~ 2 — ap 2 L 3 ). Hence the de- 
nominator eventually becomes negative and Temple's inequality looses its validity 
if L is large enough. 

As a way out of this dilemma we divide the big box A into cubic cells of side 
length £ that is kept fixed as L — > oo. The number of cells, L 3 /£ 3 , on the other 
hand, increases with L. The N particles are distributed among these cells, and 
we use ( |2.38D , with L replaced by £, N by the particle number, n, in a cell and 
p by n/£ 3 , to estimate the energy in each cell with Neumann conditions on the 
boundary. This boundary condition leads to lower energy than any other boundary 
condition. For each distribution of the particles we add the contributions from the 
cells, neglecting interactions across boundaries. Since v > by assumption, this 
can only lower the energy. Finally, we minimize over all possible choices of the 
particle numbers for the various cells adding up to A. The energy obtained in this 
way is a lower bound to Eq (A, L) , because we are effectively allowing discontinuous 
test functions for the quadratic form given by H^. 

In mathematical terms, the cell method leads to 

(2.40) E (N,L)/N > (pi^M^CnEoinJ) 

n>0 
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where the infimum is over all choices of coefficients c n > (relative number of cells 
containing exactly n particles), satisfying the constraints 



(2.41) 



E< 



i, 



E 



c n n 



pt 



The minimization problem for the distributions of the particles among the cells 
would be easy if we knew that the ground state energy Eo(n,£) (or a good lower 
bound to it) were convex in n. Then we could immediately conclude that it is best 
to have the particles as evenly distributed among the boxes as possible, i.e., c„ 
would be zero except for the n equal to the integer closest to p£ 3 . This would give 



(2.42) 



Eq(N,L) 

N 



> 4-rrpap (1 - £(p, £, R, e)) 



i.e., replacement of L in ( 2.37 ) by £, which is independent of L. The blow up of £ 
for L — > oo would thus be avoided. 

Since convexity of Eq(ti,£) is not known (except in the thermodynamic limit) 
we must resort to other means to show that n — 0(p£ 3 ) in all boxes. The rescue 
comes from superadditivity of E (n,£), i.e., the property 



(2.43) 



E (n + n', £) > E (n,£) + E (n', £) 



which follows immediately from v > by dropping the interactions between the n 
particles and the n' particles. The bound ( p. 43] ) implies in particular that for any 
n,p£ N with n > p 



(2.44) 



E(n,£)>[n/p]E(pJ)>—E(p,£) 

2p 



since the larg est in teger [n/p] smaller than n/p is in any case > n/(2p). 

The way ( 2.44| ) is used is as follows: Replacing L by £, N by n and p by n/£ 3 
in ( 2.37 ) we have for fixed R and e 



(2.45) 



E (n,£)>^p n (n-l)K(n,£) 



with a certain function K(n,£) determined by (2.38). We shall see that K is 
monotonously decreasing in n, so that if p £ N and n < p then 



(2.46) 



E (n,£)>^p n (n-l)K(p,£). 



We now split the sum in ( |2.4C| ) into two parts. For n < p we use (2.46), and for 
n > p we use (2.44) together with (|2.46) for n = p. The task is thus to minimize 



(2.47) 



2J c n n(n - 1) + I E CnU (P ~ 1 ) 

n<p n>p 



subject to the constraints (2.41). Putting 

(2.48) k := p£ 3 and t:=^c n n<k 



n<p 



we have J2 n >v Cn U = k — t, and since n(n — 1) is convex in n, and J2 n <p c ™ — ■"■ * ne 
expression ( [2.47 ) is 

(2.49) >t{t-l) + \{k-t){p-l). 
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We have to minimize this for 1 < t < k. If p > Ak the minimum is taken at t = k 
and is equal to k(k — 1). Altogether we have thus shown that 

(2.50) E (N,L) ^ ^^ p / _ J_^ R ^ pl ^ ^ 

What remains is to take a closer look at K(4p£ 3 ,£), which depends on the 
parameters e and R besides £, and choose the parameters in an optimal way. /From 
( 2.38)) and (|2.39) we obtain 



K(n,£) = (l-e)(l-^) 3 (l + fp(l-^(R 3 -R 3 )))- 1 



(2.51) x (1-- 



tt (R 3 - Rl){e£- 2 - ia£- 3 n(n - 1)) 



The estimate (2.45) with this K is valid as long as the denominator in the last 
factor in ( |2.51 ) is > 0, and in order to ha ve a f ormula for all n we can take as a 



trivial lower bound in other cases or when (2.45) is negative. As required for (2.46) 



K is monotonously decreasing in n. We now insert n — Ap£ 3 and obtain 

K{Ap£ 3 ,£) > (l-e)(l-^) 3 (1 + (const. )Y{£/ a) 3 (R 3 - R 3 )/ '£ 3 

£ 3 (const. )Y 



(2 ' 52) X \} {R 3 - R 3 ) (e(a/£) 2 - (const.)Y 2 (£/a) 3 ) 

with Y = Airpa 3 /3 as before. Also, the factor 

(2.53) (l - -^j = (1 - (const.)}- 1 ^) 3 ) 



in (2.5C) (which is the ratio between n(n — 1) and n ) must not be be forgotten 



We now make the ansatz 

(2.54) e~Y a , a/£~YP, (R 3 - R 3 )/£ 3 - Y~> 

with exponents a, (3 and 7 that we choose in an optimal way. The conditions to be 
met are as follows: 

• e(a/£) 2 — (const. )Y 2 (£/a) 3 > 0. This holds for all small enough Y, provided 
a + 5/3 < 2 which follows from the conditions below. 

• a > in order that e — > for Y — ► 0. 

• 3/3 - 1 > in order that Y- 1 {a/£f -> for for Y -> 0. 

• 1 - 3/3 + 7 > in order that Y{£/a) 3 {R 3 - Rl)/£ 3 -» for for Y -> 0. 

• 1 — a — 2/3 — 7>0to control the last factor in ( |2.52| ) . 

Taking 

(2.55) a = 1/17, /3 = 6/17, 7 = 3/17 
all these conditions are satisfied, and 

(2.56) a = 3/3-l = l-3/3 + 7 = l-a-2/3-7= 1/17. 

It is also clear that 2R/£ ~ Y 1 ! 3 = F 1 / 17 , up to higher order terms. This completes 
the proof of Theorems 3.1 and 3.2, for the case of potentials with finite range. By 



optimizing the proportionality constants in (2.54) one can show that C = 8.9 is pos- 



sible in Theorem 1.1 10] . The extension to potentials of infinite range decreasing 
faster than 1/r 3 at infinity is obtained by approximation by finite range potentials, 
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controlling the change of the scattering length as the cut-off is removed. See Ap- 
pendix B in [9] for details. A slower decrease than 1/r 3 implies infinite scattering 
length. □ 

The exponents ([2.55) mean in particular that 



(2.57) a < i? < p- 1 / 3 < J? < {paY 1/2 , 

whereas Dyson's method required R ~ p -1 ' 3 as already explained. The condition 
p -1 / 3 ^ i j g required in order to have many particles in each box and thus n(n— 1) ~ 
n 2 . The condition I <C (pa) -1 / 2 is necessary for a spectral gap gap ^> eo(p) in 
Temple's inequality. It is also clear that this choice of t would lead to a far too 
big energy and no bound for eo(p) if we had chosen Dirichlct instead of Neumann 
boundary conditions for the cells. But with the latter the method works! 



3. THE DILUTE BOSE GAS IN 2D 

The two-dimensional theory, in contrast, began to receive attention only rela- 
tively late. The first derivation of the asymptotic formula was, to our knowledge, 



done by Schick [19], as late as 1971! He found 

47T/i/9 



(3.1) e (p) 



\Hpa?)\- 



The scattering length a in ( |3.l| ) is defined using the zero energy scattering 
equation ( |2.3| ) but instead of ip(r) ~ 1—a/r we now impose the asymptotic condition 
ip(r) ~ ln(r/a). This is explained in the appendix to p|. 

Note that the answer could not possibly be eo(p) — 4irp,pa because that would 
be dimensionally wrong. But eo(p) must essentially be proportional to p, which 
leaves no room for an a dependence — which is ridiculous! It turns out that this 
dependence comes about in the ln(pa 2 ) factor. 

One of the intriguing facts about (pTj) is that the energy for N particles is not 
equal to N(N — l)/2 times the energy for two particles in the low density limit — 
as is the case in three dimensions. The latter quantity, Eq(2, L), is, asymptotically, 
for large L, equal to 87rpL~ 2 [m(£ 2 /a 2 )]~\ Thus, if the N(N - l)/2 rule were 



to apply in 2D, (3.1) would have to be replaced by the much smaller quantity 
Airpip [ln(L 2 /a 2 )] . In other words, L, which tends to oo in the thermodynamic 
limit, has to be replaced by the mean particle separation, p~ x / 2 in the logarithmic 
factor. Various poetic formulations of this curious fact have been given, but it 
remains true that the non-linearity is something that does not occur in more than 
two-dimensions and its precise nature is hardly obvious, physically. This anomaly 
is the main reason that the 2D result is not a trivial extension of Q . 

The proof of ( |3. 1[ ) is in S. The (relative) error terms to (3.1) givein in are 
| ln(pa 2 )| _1 for the upper bound and | ln(pa 2 )|~ 1//5 



To prove (|3.l|) the essential new step is to modify Dyson's lemma for 2D. The 



rest of the proof parallels that for 3D. The 2D version of Lemma 2.5 is 



Lemma 3.1. Let v(r) > with finite range Rq. Let U(r) > be any func- 
tion satisfying J U (r) ln(r / a)rdr < 1 and U(r) = for r < Rq. Let B C R 3 be 
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star shaped with respect to (e.g. convex with €E B). Then for all differentiable 
functions ip 

(3.2) f [ M |VV;| 2 + £t#| 2 ] > I* I U\i>\ 2 . 

jb Jb 



4. BOSE-EINSTEIN CONDENSATION 

Let us comment very briefly on the notion of Bose-Einstein condensation (BEC) . 
Given the normalized ground state wave function ^o(xi, • ■ ■ , £/v) we can form the 
one-body density matrix which is an operator on L 2 (R") (n = 2 or 3) given by the 
kernel 

7(f , y) = N I %o(x, &2, ■ ■ ■ , XN^oiy, x 2 ,-.., x N )dx 2 ■ ■ • dx N . 

Jbox«-i 

Then f j(£, x)dx = Trace(7) = N. BEC is the assertion that this operator has 
an eigenvalue of order N . Since 7 is a positive kernel and, hopefully, translation 
invariant in the thermodynamic limit, the eigenfunction belonging to the largest 
eigenvalue must be the constant function (volume) -1 ' 2 . Therefore, another way to 
say that BEC exists is that 

j(x, y)dxdy — 0(JV). 

Unfortunately, this is something that is frequently invoked but never proved — 
except for one special case: hard core bosons on a lattice at half-filling (i.e., N — 



half the number of lattice sites). The proof is in 20|. 

The problem remains open after about 70 years. It is not at all clear that BEC 
is essential for superfluidity, as frequently claimed. Our construction in section |2j 
shows that BEC exists on a length scale of order p^ 1 / 3 }^ -1 / 17 which, unfortunately, 
is not a 'thermodynamic' length like volume 1 ' 3 . 

5. GROSS-PITAEVSKII EQUATION FOR TRAPPED BOSONS 

In the recent experiments on Bose condensation, the particles have to be con- 
fined in a cold 'trap' instead of a 'box' and we are certainly not at the 'thermo- 
dynamic limit'. For a 'trap' we add a slowly varying confining potential V, with 
V(x) — > 00 as \x\ — ► 00 . The Hamiltonian becomes 

N 

(5.1) H = J2 -H&4 + V(xi) + Yl v (\Zi-*j\) 

i=\ l<i<j<N 

If v = 0, then 

*oO?i,...,£jv) =riiLi$o(&) 
with $o = normalized ground state of —A + V(x) with eigenvalue = A. 

The idea is to use the information about the thermodynamic limiting energy 
of the dilute Bose gas in a box to find the ground state energy of ( |5.l| ) . This has 
been done with Robert Seiringer and Jakob Yngvason |3|, Q . 

As we saw in Sections @ and there is a difference in the p dependence between 
two and three dimensions, so we can expect a related difference now. We discuss 
3D first. 
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5.1. THREE-DIMENSIONS. Associated with the quantum mechanical ground 
state energy problem is the GP energy functional 

(5.2) £ GP [$]=/ ( M |V$| 2 + y|$| 2 +4^a|$| 4 )d 3 x 

with the subsidiary condition 

|$| 2 = A 



lR 3 

and corresponding energy 

(5.3) E GP (N,a)=mi{£ GP [^} : / |$| 2 = TV} = £ GP [$ GP ]. 

As before, a is the scattering length of v. 

It is not hard to prove that for every choice of the real number N there is a 
unique minimizer for $> for £ . 

Relation of £ GP and E : If v = then clearly $ GP = y/~N <f> , and then 
£ GP = NX = Eq. In the other extreme, if V(x) = for x inside a large box of 
volume L 3 and V(x) = oo otherwise, then we take $ GP « yj N/ L 3 and we get 
E GP (N) — 4irp,aN 2 /L 3 — previous, homogeneous Eq in the low density regime. 
(In this case, the gradient term in £ GP [$] plays no role.) 

In general, we expect that E GP = Eq in a suitable limit. This limit has to be 
chosen so that all three terms in £ GP [$] make a contribution. It turns out that 
fixing Na is the right thing to do (and this can be quite good experimentally since 
Na can range from about 1 to 1000). 

Theorem 5.1. If ai = Na is fixed, 

r Eo(N,a) 
N™ocE GP (N,a) 

Moreover, the GP density is a limit of the QM density (convergence in weak 
Li-sense): 

lim ^pg(*)H*&.(2)| a • 

N—>oo 1\ 

We could imagine, instead, an N —> oo limit in which a 3> N^ 1 , i.e. a\ — » oo, 
but still pa 3 — > 0, where p is the average density, given by 

(5-4) p=lf(<P GP )\ 



N . 

In this case we simply omit the J |V$| 2 term in £ GP ; this theory is usually called 
'Thomas-Fermi theory', but it has nothing to do with the fermionic theory invented 
by Thomas and Fermi in 1927. It is appropriate for the case in which Na is much 
bigger than 1, e.g., Na w 1000. 

Proof. Outline: Getting an upper bound for E in terms of E GP is relatively 
easy, as before. The problem is the lower bound. 

One might suppose that one decomposes R 3 into small boxes as before, with 
Neumann b.c. on each box, and in each box one approximates V(x) by a constant. 

This will NOT work, even if v — 0, because all the particles will then want to 
be in the box with the smallest value of V . The gradient term will vanish and we 
will not get E GP = NX. 
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The trick is to write the quantum \&o as 

*o=niL 1 * GP (^)^i > .--»*v). 

This leads to a variational problem for F instead of ^o- Partial integration and 
the variational equation for <I> GP lead to the replacements: 

External potential: V(x) -> -87ra|$ GP (£)| 2 

and 

Measure: \[f =x dxi -» Y\f =1 \^ GV {x^dx,. 
We have to show that, with these replacements, the energy is bounded below by 
—4:TTfia f |<& GP | 4 (up to small errors). Now we can effectively use the Neumann box 
method on this F problem. 

There are still plenty of technical difficulties, but the back of the problem has 
been broken by recasting it in terms of F instead of "Jo- d 



5.2. TWO-DIMENSIONS. In view of the two-dimensional result (gT|) we 



might suppose that we must replace 47r/ia$ 4 in (5.2) by 47r/i<I> /| ln($ 2 a 2 )|, as 



suggested by (3.1). This is indeed correct, and a limit theorem for the energy and 



density is obtained that is similar to Theorem 5.1. There is a difference, however. 
This time one has to fix the ratio Nf In a instead of Na. 

Since the logarithm is such a slowly varying function it turns out that one can 
still get the correct limit and yet simplify the problem by replacing | ln(<f> 2 a 



2„2M-1 



the functional by the constant | ln(,oa 2 )| 1 , where p is the average density in (5.4) 
but for 2D, of course. 



6. THE CHARGED BOSE GAS 

The setting now changes abruptly. Instead of particles interacting with a short- 
range potential v(Si — Xj) they interact via the Coulomb potential 

v(\xi - Xj\) = \xi - Xj| _1 

(in 3 dimensions). There are N particles in a large box of volume L 3 as before, 



with p = N/L 3 . We know from the work in [21| that a nice thermodynamic limit 
exists for eo = Eq/N. 

To offset the huge Coulomb repulsion (which would drive the particles to the 
walls of the box) we add a uniform negative background of precisely the same 
charge, namely density p. Our Hamiltonian is thus 

iV 

(6.1) H = -nJ2Ai-V(x*i)+ Y, v(Si-Xj) + C 

z=l l<i<j<N 

with 

V{x) =p I \x-y\- l d 3 y and C ' = \p f V{x)d 3 x . 

Despite the fact that the Coulomb potential is positive definite, each particle 
interacts only with others and not with itself. Thus, Eq can be (and is) negative 
(just take W =const). This time, large p is the 'weakly interacting' regime. 

Another way in which this problem is different from the previous one is that 
perturbation theory is correct to leading order. If one computes ('J, H^f) with 
^ =const, one gets the right first order answer, namely 0. It is the next order 
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in 1/p that is interesting, and this is entirely due to correlations. In 1961 Foldy 
]22f calculated this correlation energy according to the prescription of Bogolubov's 
1947 theory. That theory was not exact for the dilute Bose gas, as we have seen, 
even to first order . We are now looking at second order, which should be even 
worse. Nevertheless, there was good physical intuition that this calculation should 
be asymptotically exact. It is! and this was recently proved with Jan Philip 
Solovej §. 

The Bogolubov theory states that the main contribution to the energy comes 
from pairing of particles into momenta k,—k and is the bosonic analogue of the 
BCS theory of superconductivity which came a decade later. I.e., \&o is a sum of 
products of germs of the form exp{ik ■ (xi — Xj)}. 

Foldy's energy, based on Boglubov's ansatz, has now been proved. His calcula- 
tion yields an upper bound. The lower bound is the hard part, and Solovej and I do 
this using the decomposition into 'Neumann boxes' as above. But unlike the short 
range case, many complicated gymnastics are needed to control the long range |a?| _1 
Coulomb potential. The two bounds agree to leading order in p, namely p 1 ' (al- 
though there is an unimportant technical point that different boundary conditions 
are used for the two bounds). 

Theorem 6.1. For large p 
(6.2) eo(p) ^0.402 r s ~ 3 / 4 ^ 

where r s = (3/47rp) 1 / 3 e 2 m/?i 2 . 

This is the first example (in more than 1 dimension) in which Bogolubov's 
pairing theory has been rigorously validated. It has to be emphasized, however, 
that Foldy and Bogolubov rely on the existence of Bose-Einstein condensation. We 
neither make such a hypothesis nor does our result for the energy imply the existence 
of such condensation. As we said earlier, it is sufficient to prove condensation in 
small boxes of fixed size. 

Incidentally, the one-dimensional example for which Bogolubov's theory is 
asymptotically exact to the first two orders (high density) is the repulsive delta- 



function Bose gas |17| 



To appreciate the —p' nature of (|6.2|), it is useful to compare it with what 



one would get if the bosons had infinite mass, i.e., the first term in (6.1) is dropped 



Then the energy would be proportional to — p 1 / 3 as shown in |21|. Thus, the effect 



of quantum mechanics is to lower | to j . 

It is supposedly true that there is a critical mass above which the ground state 
should show crystalline ordering (Wigner crystal), but this has never been proved 
and it remains an intriguing open problem, even for the infinite mass case. Since 
the relevant parameter is r s , large mass is the same as small p, and is outside the 
region where a Bogolubov approximation can be expected to hold. 

Another important remark about the — p 1 / 4 law is its relation to the — N 7 / b law 
for a iwo-component charged Bose gas. Dyson |23J proved that the ground state 
energy for such a gas was at least as negative as — (const)TV 7 ' 5 as N — * oo. Thus, 
thermodynamic stability (i.e., a linear lower bound) fails for this gas. Years later, 
a lower bound of this —N 7 ^ 5 form was finally established in [24], thereby proving 



that this law is correct. The connection of this — iV 7//5 law with the jellium — p 1 ^ 4 



law (for which a corresponding lower bound was also given in ]24[]) was pointed out 
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by Dyson [23] in the following way. Assuming the correctness of the — p 1 / 4 law, one 
can treat the 2-component gas by treating each component as a background for the 
other. What should the density be? If the gas has a radius L and if it has N bosons 
then p = NL~ 3 . However, the extra kinetic energy needed to compress the gas to 
this radius is NL~ 2 . The total energy is then NL~ 2 — Np 1 / 4 , and minimizing this 
with respect to L leads to the —N 7 / 5 law. A proof going in the other direction is 
in & 



A problem somewhat related to bosonic jellium is fermionic jellium. Graf and 



Solovej 1 25 1 have proved that the first two terms are what one would expect, namely 

e{p) = C TFP ^-C D p^ + o{p i /% 

where Ctf is the usual Thomas-Fermi constant and Cd is the usual Dirac exchange 
constant. 

Let us conclude with a few more details about some of the technicalities. 

6.1. FOLDY'S CALCULATION, PAIRING THEORY AND IDEAS 
IN THE RIGOROUS PROOF. Foldy uses periodic boundary conditions for 
—A, so the problem is on a torus. In order to make the Coulomb potential periodic 
and to take care of the background he replaces \x — y\~ l by 

^L- 3 H- 2 exp(zp(x-y)), 

p^O 

where the sum is over 'periodic momenta', p. (Note the p ^ condition, so the 
spatial average of this potential is 0). 
Foldy's Hamiltonian is 

N 

H' = J2 -\^ +Y.J2 i_3 br 2 eMiP&i - Xj), 

i=l i<j p^O 

in which the p ^ condition is supposed to make up for the background that is 
not explicitly included in H' . It is 'physically clear' that this device works, but a 
rigorous proof of it is not easy. 

The next step is to use the second quantization formalism, which is the one 
used by Bogolubov and which is a very convenient bookkeeping device (but it has 
to be noted that it is no more than a convenient device and it does not introduce 
any new physics or mathematics). 

H' = ^2 \p\ 2 a*p a P + Yl L ~ 3 \P\~ 2 5Z a t+p a *q-p a q a k, 
p p^O k,q 

where the a p operators satisfy the usual bosonic commutation relations. An im- 
portant observation is that since p ^ there are no terms with 3 or 4 oj. This is 
different from the situation with the usual short range potential treated in Section 
H and it means that the leading term in perturbation theory vanishes. Everything 
now comes from the terms with two a . 

6.1.1. The Bogolubov approximation: The motivation is Bose condensation: Al- 
most all particles are in the state of momentum p = created by a^. Thus: 
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Step 1 in the Bogolubov approximation: Keep only quartic terms with precisely 
two a (and ignore terms with one or no a ). We obtain the reduced Hamiltonian 



p p#o 



a p a a p ao 



+a* a*_ p a a- p + a p a*_ p aoao + aQCiQa p a- p 



Step 2 in the Bogolubov approximation: Replace the operators a by the number 
N. We then obtain: 



H '" = Yl \p\ 2a * P a p+p\p\~ 

p#0 



a p a p + a_ p a- p 



+a p a_ p + a p a- p 

This is a quadratic Hamiltonian and can be diagonalized by completing the square: 

H'" = 2_, Ap(a P + P P a- p ){a p + f3 p a*_ p ) 
p 
+A p (a*_ p + f3 p a p )(a- p + f3 p a* p ) 



-2 £4,/?; 



p#0 



The last term comes from the commutator [a p , a*] = S pq . 

A P {l + P 2 p ) = \\P\ 2 + P\P\- 2 
2A p f3 p = p\p\- 2 

The ground state energy is given by the last term above. 

e = ^lim ~ J2 A A = I A pPp = Cp P 5/i - 

In this approximation the ground state wave function tp satisfies 

(a p + l3 p aL p )il) = 0, 

for all p^O. 

In the original language (in which oq an operator) this corresponds to a function 
of the form 



v> 



l + ^2f{Xi ~Xj) 



l<] 



+c ^ f(xi-Xj)f(xi- Xk) + ■■■ 

different 

where /(p) = (3 P . In fact, f(p) = G(|p| 4 /p), with G independent of p. Thus / varies 
on a length scale p -1 ' 4 (which is the typical inter pair distance). 
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6.1.2. Ideas in the rigorous proof: As in the short range case in Section 0, there 
is no need to prove Bose condensation globally. It is enough to do so on a short 
scale. 

• Localize particles, by means of Neumann bracketing, in "small" boxes of 
size £. The constant function (i.e., the 'condensate') is not affected by this 
localization since the constant function in a small box satisfies Neumann 
boundary conditions. The function / discussed above is not affected very 
much if £ 3> p~ x l 4 . We choose £ close to p~ x l 4 . 

• We control the Coulomb interaction between boxes by using an averaging 



method in [24|. The error in neglecting intercell interactions can be shown 
to be dominated by N/£ < Np 1 / 4 . 

We establish condensation on the scale £ by noting that the first non-zero 
Neumann eigenvalue is ~ £~ 2 . If N + is the expected number N + of particles 
not in the condensate in the "small box" , their energy is bounded from 
below by N + £~ 2 ~ N + p 1 / 2 . If this is to be consistent with the total energy 
—Np 1 / 4 we should expect to have N + <C Np~ 1 / 4 , i.e., local condensation 
exists. 

This is established by means of a bootstrap procedure. 
Having established local condensation one starts the hard work of proving 
that the Bogolubov 'quadratic' approximation really gives the leading term 
in the energy. There is however a new difficulty that arises from the fi- 
nite size of the small boxes in which we are working. Neumann, and not 
periodic boundary conditions must be used. Since we no longer have 'mo- 
mentum conservation' in the small boxes, diagonalization of the Hamiltonian 
by 'completing the square' is not the simple algebraic problem it was before. 
Nevertheless, it can all be carried to a succesful conclusion 
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